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CO ' Abstract 

We report on results of astrometric observations of H2O masers in the "water fountain" source 
IRAS 18286-0959 (118286) with the VLBI Exploration of Radio Astrometry (VERA). These observations 
£/3 ' yielded an annual parallax of IRAS 18286—0959, ir = 0.277 ±0.041 mas, corresponding to a heliocentric 

distance of D = 3.61±£ !g kpc. The maser feature, whose annual parallax was measured, showed the ab- 
solute proper motion of (fi a ,^s) = (— 3.2 ± 0.3, — 7.2 ± 0.2)[mas yr -1 ]. The intrinsic motion of the maser 
feature in the internal motions of the cluster of features in 118286 does not seem to trace the motion of the 
bipolar jet of 118286. Taking into account this intrinsic motion, the derived motion of the maser feature is 
roughly equal to that of the maser source 118286 itself. The proximity of 118286 to the Galactic midplanc 
(z sslO pc) suggests that the parental star of the water fountain source in 118286 should be intermediate- 
mass AGB/post-AGB star, but the origin of a large deviation of the systemic source motion from that 
expected from the Galactic rotation curve is still unclear. 
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• 1. Introduction AGB stars (e.g., Imai 2007). 

' Here we report on the measurement of an annual par- 
Energetic mass loss (rate up to lO _3 M0yr _1 ) has been allax of H2O masers in IRAS 18286—0959 (hereafter ab- 
observed from dying stars such as asymptotic giant branch breviated as 118286) with the VLBI Exploration of Radio 
(AGB) and post-AGB stars. The spatio-kinematics of Astrometry (VERA). The VERA astrometry for 118286 
the mass-loss flows have been investigated in great de- was conducted in one of the VERA key science projects, 
tail by observations of maser emission such as SiO, H2O, which focuses on H2O maser sources at intermediate dis- 
. ' an( l OH using very long baseline interferometry (VLBI). tances (2-5 kpc from the Sun) in order to extend the dis- 
Thcse maser sources are associated with circumstellar en- tance scale for the VERA parallax measurements. The 
$H • velopes (CSEs) of AGB stars, but in rare cases with highly result of the 118286 astrometry could be used as one of 
collimated bipolar jets, so called "water fountains," and the data points for exploring the Galactic dynamics if the 
CSE remnants around post-AGB stars or central objects motion of 118286 closely follows that of the Galactic ro- 
of planetary nebulae (PNe). In the latest decade, high tation (i.e. Reid et al. 2009; Honma et al. 2012). 118286 
precision astrometry in very long baseline intcrfcromery is a member of the class of water fountains, whose spatio- 
(VLBI) has enabled the measurement of trigonometric kinematics of H2O masers have been investigated in de- 
parallax distances and proper motions of maser sources, tail. Yung et al. (2011) suggest that most of the H2O 
Even in the rare cases mentioned above, the heliocentric masers in 118286 are associated with "double helices" of 
distances and three-dimensional secular motions of the highly collimated bipolar jets with speeds of ~100 km s _1 . 
water fountains, pre-PNe, and PNe have been measured, The formation of the double helices is expected from dis- 
in those cases where H2O masers could be detected (Imai continuous mass ejection from a moving star on a pe- 
et al. 2007c; Imai et al. 2011b; Tafoya et al. 2011). These riod of ^30 yr. The annual parallax measurement for 
kinematic approaches have contributed to estimating the 118286 is the second case of parallax measurements for 
physical parameters of these parental stars and their evo- the water fountains after that for IRAS 19134+2131 (Imai 
lutionary properties. Taking into account their locations et al. 2007c). This paper discusses the galactic kincmatical 
and 3D motions in the Milky Way, it has been suggested properties of 118286 in Sect. 3.2. However, the existence 
that these stars should be intermediate- mass AGB/post- of a low- velocity "equatorial flow" should be taken into 
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account when one interprets the observed secular motion 
of the maser feature in the cluster of features in 118286 
(Imai 2007). The argument of the equatorial flow will be 
described in a separate paper in more detail. 

2. Observations and data reduction 

The VERA observations of the 118286 H 2 ( J k ^k + = 
612 — > 523 , 22.2 3 5 80 GHz) masers were conducted at 16 
epochs from 2007 October to 2009 September. Table 1 
gives a summary of these observations, maser source map- 
ping, and astrometry. Out of 16 epochs, 14 produced 
scientifically meaningful output. At each epoch, the ob- 
servation was made for 6.7 hr in total. 118286 was ob- 
served together with the fringe-phase and position refer- 
ence source, ICRF J183220.8-103511 (hereafter abbrevi- 
ated as J1832), separated by 0°.67 from 118286, simul- 
taneously using VERA's dual-beam system. They were 
tracked for 30-35 min out of every 40 min, followed by 
scans on other band-pass calibrator sources. The received 
signals were digitized in four quantization levels, then di- 
vided into 16 base band channels (BBCs) with a band- 
width of 16 MHz each. One of the BBCs collected data 
from 118286 at the frequency band including the H2O 
maser emission, while other BBCs from J1832 spanning 
a total frequency band range of 480 MHz. The BBC out- 
puts had a recording data rate of 1024 Mbits s _1 . The 
data correlation was made with the Mitaka FX correla- 
tor with a correlation accumulation period of 1 s. The 
correlation outputs consisted of 512 and 32 spectral chan- 
nels for the H2O maser and reference continuum emission, 
respectively. The former corresponds to a velocity spac- 
ing of 0.42 km s" 1 , which is narrow enough to resolve an 
H2O maser feature (corresponding to a maser gas clump), 
which consists of two or more spectral channel components 
called maser spots. 

Data reduction was mainly made with the National 
Radio Astronomy Observatory (NRAO) Astronomical 
Image Processing System (AIPS) package. For astrom- 
etry, we need special procedures described as follows 
(see also e.g., Honma et al. 2007; Imai et al. 2007b). 
Firstly, delay-tracking was repeated for the correlated 
data using better delay-tracking models calculated with 
the software equivalent to the CALC9 package developed 
by the Goddard Space Flight Center/NASA VLBI group. 
Throughout the whole data analysis, we adopted the coor- 
dinates of the delay-tracking center: aj2ooo 

= 18 h 31 m 22 s 

■934, ^ooo = -09°57'21".70 for 118286 and 

a/2000 =18 h 32 m 20 s .836521, <5 J20 oo = -10°35'11". 20055 
for the reference source J1832. The delay-tracking 
solutions include residual delay contributions from the 
atmosphere, which were estimated using the global 
positioning system (GPS) data (Honma et al. 2008b). 
Secondly, differences in instrumental delays between two 
signal paths in the dual beam system were calibrated 
using the differential delays, which were measured using 
artificial noise signals injected into the two receivers at 
the same time. The measurement accuracy has been 



improved for observations since 2007 by installing one 
to four artificial noise sources (Honma et al. 2008a). 
Thirdly, fringe-fitting and self-calibration were performed 
using the continuum source data, whose solutions were 
applied in the data analysis of maser emission. Only the 
solutions in the BBC at the same frequency as that for 
the H2O maser emission were valid for the fringe-phase 
calibration. The accuracy of the phase-fluctuation 
compensation and the coherence of the integration is 
dependent on weather condition, particularly the humid- 
ity, which is seasonally variable. These effects affect the 
final astrometric accuracy. Finally, image cubes of the 
maser source were obtained by deconvolution through 
the CLEAN algorithm with a typical synthesized beam 
of 0.9x2.5 in milliarcseconds (mas) in the case of full 
operation of VERA's four antennas. Because the CLEAN 
deconvolution was performed by automatically selecting 
local peaks of brightness in each of channel maps without 
smaller CLEAN boxes, side lobes of the synthesized 
beam pattern affected the final CLEAN image cube 1 . 
Each maser spot (or velocity component) was identified 
as a Gaussian brightness component using the AIPS task 
SAD. 

In order to obtain higher quality maser maps, we also 
adopted a normal procedures for maser source mapping. 
Fringe-fitting and self-calibration were performed using 
a spectral channel that contain bright maser emission. 
Column 4 in table 1 gives the local-standard-of-rest (LSR) 
velocity of the spectral channel selected as phase- and 
position-reference. The obtained solutions of calibration 
were applied to the data in all spectral channels. 

3. Results 

3.1. The distribution and proper motions of H2O maser 
features in IRAS 18286-0959 

Column 6 of table 1 gives the numbers of maser features 
identified in the individual epochs. Figure 1 shows the dis- 
tributions of H2O maser features at six out of 14 epochs. 
The maser feature distribution was highly variable, with 
only a small fraction of the maser features having lifetimes 
longer than half year. Note that the VLBA observations 
conducted by Yung et al. (2011) and the VERA observa- 
tions reported in this paper, respectively, lasted for simi- 
lar seasons, 2008 April-2009 May and 2007 October-2009 
September, respectively. The latter observations identi- 
fied most of the same maser features as found in Yung 
et al. (2011). However, there exists some discrepancy of 
maser feature detections between these observations be- 
cause of the different image sensitivities 2 and, more im- 

1 A tool for image deconvolution with automatic CLEAN 
box setting is currently in construction and will be 
used in future image synthesis. This tool works in a 
Python/ParselTonguc script. ParselTonguc is a library 
developed in JIVE and can control AIPS tasks from the 
outside of AIPS POPS using Python scripts. See this URL. 
http: / /www .jive. nl/dokuwiki/doku.php?id=parscltonguc:parscltongue 

2 The H2O maser features of 118286 were resolved out in longer 
baselines of the VLBA as described in table 1 of Yung 
et al. (2011). Also the total integration time of 118286 was 
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portantly, different epochs of observations. Although our 
VERA observations could not trace the maser feature that 
was selected as position-reference in the VLBA observa- 
tions (at an LSR velocity Vlsr— 52 km s _1 in the eastern 
side of the maser distribution) over a full year, we newly 
identify another maser feature that survived over one year 
in the western side of the maser distribution. This maser 
feature (IRAS 18286-0959:12013- U, see table 2), has a 
similar LSR velocity and is denoted by a black opened 
circle in figure 1. The maser spots included in this maser 
feature was used for the annual parallax measurement as 
described in subsection 3.2. 

The relative proper motions of maser features with 
respect to the position-reference feature were identi- 
fied for the features detected at three or more epochs. 
Some proper motions were also removed from considera- 
tion because their LSR velocity drifts were higher than 
10 km s _1 yr _1 . Thus only 50 maser proper motions were 
measured although there existed a few hundred identified 
features. Because there was no single maser feature that 
was detected at all epochs, different maser features were 
selected as position reference in the first nine and final 
five epochs. Table 2 gives the parameters of the maser 
features whose proper motions were measured. Figure 2 
shows the spatial distribution of the relative maser proper 
motions. The spatio-kinematical structure roughly repro- 
duces a fast bipolar outflow found by Yung et al. (2011), 
but it is difficult to recognize the double helix jets as seen 
by Yung et al. (2011). The VLBA identification of maser 
proper motions seems to have biases to high velocity mo- 
tions. It is noteworthy that, in the present result, there 
exist a group of maser features moving in slow veloci- 
ties (V ^30 km s _1 ) near the center of the maser fea- 
ture distribution. They were moving perpendicularly to 
the high velocity (V ^100 km s _1 ) motions of maser fea- 
tures in the north-south direction. This implies the exis- 
tence of an "equatorial flow" (Imai 2007), different from 
the high velocity flows (V ^100 km s _1 ) that are associ- 
ated with highly collimated bipolar jets described by Yung 
et al. (2011). This equatorial flow will be described in a 
separate paper. 

Nevertheless it is still necessary to estimate the motion 
of the star itself in 118286 (or the systemic motion) with 
respect to the maser proper motion in order to discuss 
the systemic motion of 118286 in the Milky Way (see sub- 
section 4). The maser motions shown in figure 2 are well 
approximated to an outflow with a point-symmetric veloc- 
ity field. We performed the least-squares fitting method 
for the maser motion data using a radially expanding flow 
model. This method has been repeatedly described in our 
previous papers (e.g. Imai et al. 2011a; ?). Note that the 
model fitting assumes independent radial-expansion veloc- 
ities of the maser features rather than a common radial 
velocity field as a function of distance from the dyamical 
center of the flow. This is, however, sufficient for estimat- 

shortcr in the VLBA observations than that of the VERA ob- 
servations. Therefore, the difference of the array sensitivity was 
not so significant as expected from full operation of VLBA's 10 
antennas. 



ing the systemic motion of the outflow origin as one of the 
free parameters. Table 3 gives the derived fitting param- 
eters. Figure 2a displays the measured maser motion vec- 
tors after subtracting the derived systemic motion vector. 
The model-fitting result suggests that position-reference 
feature (IRAS 18286-0959:12013-14) itself may belong to 
the equatorial flow mentioned above. 

3.2. The annual parallax distance to IRAS 18286—0959 

Column 7 of table 1 shows validity of astrometry for 
the annual parallax measurement in the individual epochs 
(Y/N/S). As described in subsection 3.1 and table 1, 
the maser feature IRAS 18286-0959:12013-^ was iden- 
tified in the observations between 2007 October and 2009 
March in the maser image cubes synthesized through self- 
calibration and/or phase- referencing procedures and two 
maser spots in this feature were used for the parallax mea- 
surement. Figure 3 shows the motions of these maser 
spots. The spot motions can be fitted to combination 
of an annual parallax and a constant-velocity motion as 
described later. The uncertainty of spot position pre- 
sented in figure 3 indicates only the contribution from 
thermal noise, which is underestimated if the spot is ex- 
tended. We note that Yung et al. (2011) and the data of 
Imai (2007) suggest that the H2O maser spots in 118286 
were extended. Thus, deviations of the data points from 
the modeled motions (up to 0.5 mas) are larger than the 
position errors expected from combination of contribu- 
tions from thermal noise, instrumental, and atmospheric 
phase-delay residuals (0.1 mas level, Honma et al. 2010, 
see also general formulation in Asaki et al. 2007). They 
should be mainly attributed to the variation of maser 
spot structures. Note that the temporal variation of 
maser spot/feature structure is not necessarily random 
as assumed in Honma et al. (2010) if the feature is as- 
sociated with some physical feature (e.g. a shock front, 
Imai et al. 2002) changing on a specific time scale. 
Nevertheless, random variation of a spot position around 
its ballistic motion as shown in figure 3 is expected when 
the maser is highly spatially resolved with a VLBI synthe- 
sized beam whose shape may change from one epoch to 
another as small intrinsic variation of the maser structure 
will be highly enhanced (c.f. Imai et al. 2007b). 

We attempted the least-square method for the spot mo- 
tions to fit the modeled motions each of which is composed 
of an annual parallax, a constant secular motion, and a 
position offset at the reference epoch (J2000.0) 3 . In order 
to obtain the common annual parallax of these spots, it- 
erative procedures were adopted, which is similar to the 
approach by Sanna et al. (2012). Firstly, the model fit- 
ting was performed independently for the individual spot 
motions (independent fitting). Secondly, a mean annual 
parallactic motion was subtracted from the spot motions 
to estimate only the position offsets and secular motions 
of spots (proper motion fitting). Thirdly, the systemic 
motions estimated from the derived parameters were sub- 

3 In the model fitting, all spot positions are at first converted to 
the relative positions with respect to that at the first observation 
epoch for conserving the fitting accuracy. 
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tracted from the original spot motions. Finally, the po- 
sition residuals were used to estimate a common parallax 
and position offset and proper motion residuals (combined 
fitting). The proper motion fitting was again performed 
using the estimated common parallax. Then the proce- 
dures from the proper motion to combined fitting were it- 
erated until the estimated parameters seemed converged. 

Table 4 gives the derived parameters after the final it- 
eration. Although the two maser spots were associated 
with the same maser feature and their motions are not 
completely independent, the data combination may mit- 
igate the error contribution from statistical (large) er- 
rors expected from the signal-to-noise ratios of maser 
spot detection and the random variation of spot struc- 
tures. In the model fitting, weighting with position accu- 
racy was adopted as performed in our previous analysis 
(Imai et al. 2011b). Thus we obtained an annual paral- 
lax of 118286, 7r = 0.277 ±0.041 mas, corresponding to a 
distance value of D = 3.61^0 47 kpc. The kinematic dis- 
tance to 118286 of 3.1 kpc (Dcguchi et al. 2007) is roughly 
consistent with the annual parallax distance. 

The absolute proper motion of the maser feature 
IRAS 18286-0959:12013-1^ was derived to be (fj, a ,ti s ) = 
(—3.2 ± 0.3, —7.2 ± 0.2) [mas yr -1 ] from the mean motions 
of the two spots, whose parameters are given in table 4. 
Here wc estimate the location and the three-dimensional 
motion of the 118286 system in the Milky Way. Table 
5 gives the derived parameters. 118286 is located very 
close to the Galactic midplane (z ~7 pc). On the other 
hand, the systemic secular motion of 118286 has some un- 
certainty, described as follows, which is dependent on the 
relative motion of the feature IRAS 18286-0959:12013-1^ 
with respect to the central star. If the motion of this 
maser feature is the same as that of the star (Case 1 in 
table 5), the systemic motion has a very large deviation 
(up to 100 km s _1 ) from that expected from the Galactic 
rotation curve. This is a similar case to that in which one 
adopts the relative motion of the feature estimated from 
the model fitting as described in subsection 3.1 (Case 2). 
The systemic motion only follows the Galactic rotation 
curve closely if one assumes a relative proper motion of 
the feature (—3,-5) [mas yr" 1 ], or (—51, —85) [km s _1 ] 
with respect to the star, as indicated with a magenta ar- 
row in figure 2 (Case 3). 

4. Discussion 

Through the present work, we can obtain some con- 
straint on the properties of the central star in the water 
fountain 118286. The annual parallax distance to 118286 
gives the luminosity value of 118286, ~ 1.2 x 10 4 L e , a 
factor of 1.4 higher than the value previously derived 
(8.7 x 10 3 L Q , Deguchi et al. 2007). This suggests that 
118286 has a luminosity higher than the typical value of 
AGB stars (~6 000 L©), supporting the hypothesis of 
118286 to be a higher mass AGB/post-AGB star (e.g., 
Imai 2007; Imai ct al. 2007c). Note that the high mass 
population of evolved stars such as red supcrgiants follows 
closely the Galactic rotation within 30 km s _1 (e.g. Rcid 



ct al. 2009; Asaki et al. 2010; ?) while lower mass evolved 
stars may have the large kinematical deviations due to the 
dynamical relaxation in the Milky Way during their long 
lifetime. In the case of 118286, the large deviation of its 
motion from the Galactic rotation favors the hypothesis 
of an AGB star over that of a red supergiant. However, its 
close proximitiy to the Galactic mid-plane does not rule 
out the possibility that the large kinematical deviation is 
attributed to a binary motion in 118286. SiO maser or 
millimeter continuum emission can pinpoint the location 
of the star (Imai et al. 2005; Imai 2007), but has not yet 
been detected toward 118286. Detection of a circumstcl- 
lar envelope in molecular line emission with a new facil- 
ity such as the Atacama Large Millimeter/submillimeter 
Array (ALMA) will shed light on the property of of the 
central star. In practice, it is difficult to distinguish 
the intrinsic emission of 118286 from interstellar fore- 
ground/background emission (Imai et al. 2009) although 
this object is clearly identified in the Spitzer/ GLIMPSE 
image (Deguchi et al. 2007). In addition, more intensive 
monitoring VLBI observations of the 118286 H2O masers 
are necessary in order to trace a greater number of proper 
motions of maser features including shorter-lived (a few 
months) features, enabling to elucidate the maser spatio- 
kincmatics in more detail, including a possible equatorial 
flow or an AGB envelope. 
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Table 1. Parameters of the VERA observations 



Observation 


VERA 






Beam 


i 










Astromctry 


code 


epoch 


telescopes 1 * 


(km s _1 ) 


Noisc c 


[mas 






valid?' 


/pc 


,k A 


Note 


r07296a 


2007 October 23 


MROS 


48.2 


42 


4.34x0.98, 


-47". 2 


68 


Y 


1 


5 


Calibration solutions invalid (O) 


r07328a 


2007 November 26 


MROS 


48.2 


-12 


2.35x0.91, 


-41° 1 


71 


Y 


1 


8 


r07357b 


2007 December 26 


MROS 


48.0 


31 


2.64x0.83, 


-37° 1 


(iO 


Y 


1 


3 




r08048b 


2008 February 17 


MROS 


46.7 


24 


1.63x0.88, 


-30°.2 


105 


Y 


2 


5 




r08098a 


2008 April 9 


MROS 


70.1 


37 


1.65x0.96, 


-32°5 


73 


Y 


2 


4 




r08132a 


2008 May 13 


MROS 


48.3 


40 


4.09x0.90, 


-37°8 


03 


S 


2 


3 


See comment * 


r08166a 


2008 Juno 16 


(R)(0)S 












N 






Earthquake (M), heavy rain (R) 


r08174a 


2008 June 24 


MOS 












N 






High noise temperature (O) 


r08189a 


2008 July 9 


MROS 


86.3 


90 


1.31x1.24, 


+65°6 


01 


Y 


2 


8 


Too faint detected reference source 


r08223a 


2008 August 12 


MROS 


48.1 


07 


1.59x1.18, 


-49°0 


53 


Y 


2 






r08329a 


2008 December 26 


MROS 


83.7 


50 


2.22x1.05, 


+89°.6 


97 


Y 


1 







r09019b 


2009 January 19 


MROS 


100.7 


58 


4.06x0.77, 


-38°2 


69 


Y 


1 





See comment 


r09050b 


2009 February 19 


M(R)0(S) 


103.5 


100 


2.68x0.95, 


+ 14°. 9 


34 


N 






High noise temperatures (R, S) 


r09068b 


2009 March 9 


MROS 


94.2 


40 


2.17x0.86, 


-35°9 
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High noise temperatures (O) 



a Tclescopc whose data were valid for phase-referencing mascr imaging. M: Mizusawa, R: Iriki, O: Ogasawara, S: Ishigakijima. The station with 
parentheses had some problem during the observations and affected the annual parallax measurements. 

LSR velocity at the phase-referenced spectral channel. 
c Smallcst rms noise in the emission-free spectral channel image obtained through self-calibration image synthesis, in units of mjy beam -1 . 
d Synthcsizcd beam size resulting from natural weighted visibilities, i.e. major and minor axis lengths and position angle. 
°Number of the detected mascr features in the data obtained through self-calibration. 

Y and N: valid and invalid data point for annual parallax measurement, respectively. S: valid in a special procedure, in which the maser position 
was obtained from the relative position of another maser spot whose coordinates were measured in both of the image cubes synthesized through 
self-calibration and phase-referencing procedures. 

s Peak intensity, in units of Jy beam - , of the phase-referencing-based image of the maser feature IRAS 18286—0959:12013-.?^ that contained the 
maser spots selected for the astrometry. 

A too long fringe-fitting solution interval (3 min.) was adopted. The maser image cube obtained through self-calibration did not have better 
quality than that through phase-referencing calibration, causing some of the mascr features found in the latter image cube missing in the former 
one. 
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Table 2. Parameters of the H2O maser features identified by proper motion toward IRAS 18286—0959 



Offset Proper motion 13 Radial motion c Detection at 

(mas) (mas yr"" 1 ) (km s -1 ) individual epochs 
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a H 2 maser features detected toward IRAS 18286-0959. The feature is designated as IRAS 18286-0959:12013- N, where Ms the ordinal source 

number given in this column (12013 stands for sources found by Imai ct al. and listed in 2013). 

b Rclative value with respect to the motion of the position-reference maser feature located at the map origin. 

c Rclative value with respect to the local standard of rest. 

d Mean full velocity width of the maser feature at half intensity. 
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Table 3. Parameters of the best fit 3D spatio-kinematical model of the H2O masers in IRAS 18286—0959 



Parameter Value 

AX [mas] * 24±3 

AYh [mas] a 49±4 

Vox [km s _1 ] a 13±2 

V Y [km s" 1 ] a 21±4 

Vaz [km s" 1 ] b 65 

D [kpc] b 3.6 

y^" 4~22~ 



a Rclative value respect to the maser feature IRAS 18286 — 0959:12013--/^. 
b Assumcd value. 

c The x 2_sc i uare value was obtained in the model fitting based on the least-squares method. 



Table 4. Parameters of the fitted maser spot motion 



^LSR 


X " 




fJ, Y 


7T 


[km s _1 ] 


[mas] 


[mas] [mas yr~ 1 ] 


[mas yr~ 1 ] 


[mas] 


Independent fitting 


53.4 


26.45±1.00 


55.56±0.86 -3.37±0.12 


-7.11±0.10 


0.284±0.058 


53.0 


27.57±1.20 


56.53±0.90 -3.50±0.15 


-7.24±0.11 


0.229±0.056 


Proper motion fitting (after subtracting the annual parallax b 


53.4 


25.46±1.60 


55.46±0.98 -3.26±0.19 


-7.10±0.12 




53.0 


25.10±2.22 


56.90±0.94 -3.22±0.27 


-7.28±0.11 




Combined fitting to an annual parallax and residuals of position offset and linear motion 


Combined 


1.56±0.80 


-0.06±0.63 -0.17±0.10 


0.01±0.08 


0.277±0.041 



a Position at the epoch J2000.0 with respect to the delay-tracking center (see main text). 

b The adopted annual parallax was given in the combined fitting in the converging phase of iterations. 



Table 5. Location and 3D motion of IRAS 18286-0959 in the Milky Way estimated from the VERA astrometry 



Parameter Value 

Galactic coordinates, (I, b) [deg] a (21.80, -0.13) 

Heliocentric distance, D [kpc] a 3.61±0.63 

Systemic LSR velocity, V sys [km s — ] a . . 60±5 

Hi, [kpc] b 8.0 

6o [km s _1 ] b 220 

(U Q , V ,W O ) [km s~ 1 ] c (7.5,13.5,6.8) 

20 [pc] d 16 

/>' .„ [kpc] 4.84±0.50 

z [pc] 7±1 

Case 1: fixed in the system 

V R [km s" 1 ] 66±19 

V e [km s" 1 ] 148±25 

V z [km s -1 ] -1±20 

Case 2: moving in ( — 0.76, —1.23) [mas yr _1 ] with respect 
to the system (result of the kincmatical model fitting) 

V R [km s" 1 ] 51±17 

Ve [km s" 1 ] 169±22 

V z [km s- 1 ] -2±20 

Case 3: moving in ( — 3, —5) [mas yr~ I ] 
with respect to the system 

Vr [km s _1 ] 6±13 

Ve [km s _1 ] 228±17 

V z [km s- 1 ] -7±20 

a Input value for IRAS 18286-0959. 

b Input value for the Sun in the Milky Way. 



c Motion of the Sun with respect to the local standard of rest, cited from Francis &c Anderson (2009) (c.f, Dehnen & Binney 1998). 
d Hcight of the Sun from the Galactic mid-plane, cited from Hammerslcy et al. (1995). 
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Fig. 1. Distributions of H2O maser spots in IRAS 18286—0959. A black opened circle indicates the location of the maser feature 
(IRAS 18286-0959:12013-^4) that contains maser spots measured their annual parallax, (a) On 2007 October 23. (b) On 2008 
February 17. (c) On 2008 August 12. (d) On 2009 January 19. (c) On 2009 May 11. (f) On 2009 September 9. We note that 
some maser spots that were detected in the image cube obtained through phase-referencing calibration are missing in these plots, 
especially on 2009 January 19. Such maser spots include those detected in the phase-referencing image synthesis and measured their 
annual parallax. 
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Fig. 2. Relative proper motions of H2O maser spots in IRAS 18286—0959. The origin of coordinates is set to the position- reference 
maser feature, which was selected from different features between the epochs in 2007-2008 (IRAS 18286— 0959:12013- 14) and those 
in 2009 (IRAS 18286—0959:12013-36) . Colors of maser feature indicate LSR velocities. An arrow shows the relative proper motion 
of the maser feature. The root position of an arrow indicates the location of the maser feature at the first of the epochs when 
the feature was detected. The length and the direction of an arrow indicate the speed and direction of the maser proper motion, 
respectively, (a) Proper motions identified in 2007-2008. The annual parallax and secular motions of maser spots were measured for 
those in the position-reference feature at the map origin. The systemic motion (or the motion of star), which was estimated in the 
least-squares model fitting, is subtracted from the measured proper motions. A magenta plus symbol denotes the estimated location 
of the dynamical center of the modelled outflow. A magenta arrow indicates a proper motion (—51, —85) [km s -1 ], corresponding to 
(—3,— 5) [mas yr~ x ]. If the relative motion of the position-reference feature has this motion vector with respect to the system, the 
systemic motion well follows a motion expected from the Galactic rotation curve, (b) Same as (a) but in 2009. A mean motion of 
(20, -130) [km s" 1 ] is subtracted from the individual measured motions. 
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Year since J2000.0 



Fig. 3. Motions of the 53.4 km s" 1 and 53.0 km s _1 components of H2O masers in IRAS 18286—0959 and the kinematical models 
for these motions, (a) R.A. and decl. offsets with respect to the phase-tracking center of the 53.4 km s -1 component. A filled circle 
shows the data point observed and used for the annual parallax measurement. A solid curve shows the modeled motion including an 
annual parallax and a constant velocity proper motion. An opened circle indicates the spot position expected in the model at the 
observation epoch, (b) R.A. variation of the 53.0 km s — 1 component along time. The estimated linear proper motion is subtracted 
from the observed spot position. A solid curve shows the modeled annual parallactic motion, (c) Same as (a) but for the 53.4 km s _1 
component, (d) Same as (b) but for the 53.0 km s _1 component, (e) The result of the combined annual parallax fitting. Blue 
and red data points shows those of the 53.4 km s _1 and 53.0 km s _1 components, respectively. For clarity, data points are slightly 
shifted in the horizontal axis. 



